ADA103054 


AFWL-TR -79-50 


(D  LEVEL  ^ 


AFWL-TR- 

79-50 


ELECTROMAGNETIC  CONSIDERATIONS  OF  A 
SPATIAL  MODAL  FILTER  FOR  SUPPRESSION  OF 
NON-TEM  MODES  IN  THE  TRANSMISSION-LINE 
TYPE  OF  EMP  SIMULATORS 

D.  V.  Giri 
C.  E.  Baum 
H.  Schilling 


LuTech,  Inc 
P.O.  Box  1263 
Berkeley,  CA  94701 

June  1980 


Final  Report 


s 


DTIC 

ELECTE* 
AUG  1  9  1981 J 


Approved  for  public  release;  distribution  unlimited 


AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base,  NM  87117 


3?2> 


81  8  18  085 


AFWL-TR-79-50 


This  final  report  was  prepared  by  LuTech,  Inc,  Berkeley,  California, 
under  Contract  F29601 -78-C-0080,  Job  Order  12090606  with  the  Air  Force  Weapons 
Laboratory,  Kirtland  Air  Force  Base,  New  Mexico.  Dr.  Michael  G.  Harrison  (NTYE) 
was  the  Laboratory  Project  Officer-in-Charge. 

When  US  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  a  definitely  related  Government  procurement  operation,  the 
Government  thereby  incurs  no  responsibil ity  nor  any  obligation  whatsoever,  and 
the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded 
by  implication  or  otherwise,  as  in  any  manner  licensing  the  holder  or  any  other 
person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use, 
or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  authored  by  a  contractor  of  the  United  States 
Government.  The  United  States  Government  retains  a  nonexclusive,  royalty-free 
license  to  publish  or  reproduce  the  material  contained  herein,  or  allow  others 
to  do  so,  for  the  United  States  Government  purposes. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be 
available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


ty/LJicuJ) Q] 

MICHAEL  G.  HARRISON,  PhD 
Project  Officer 


oJ.  PHILIP  CASTILLO,  PhD 
Chief,  Electromagnetics  Branch 


FOR  THE  DIRECTOR 


THOMAS  W.  CIAMBRONE 
Colonel,  USAF 

Chief,  Applied  Physics  Division 


DO  NOT  RETURN  THIS  COPY.  RETAIN  OR  DESTROY. 


UNCLASSIFIED 


SECURt-v  Z±MS\r ICATIQN  OF  this  RAGE  'Whtn  Daia  Enttrtb) 

( /^./REPORT  documentation  page 


ffi. 


AFWlVfR-79- 


1 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


;2.  GOVT  ACCESSION  NO.i 

i 


(\^03 


i.  reorien-'icatalog  numier 


X.  f  I Ti_£  and  Submit) 


ELECTROMAGNETIC  CONSIDERATIONS  OF  A  JPATIAL 
JODALilLTER  FOrSUPPRESSION  OF  NON-Tff  -MODES 
iN  THE  JRANSMISSfON^INE  TYPE  OF  .gJP  SIMULATORS, 


-ty»i  o*-RtRo»r.t  reriqo  covereo 

inaiyfeiJevt.  ) 

RERFORhuliyWnfWERORT  n 


umBER 


D.  V./Giri 
C.  E.Alaum 
H .  ^5c  hilling 


8.  CONTRACT  OR  GAAnT  number*; 

F2960  i  -7S-C-008&I  t>  *-  'J 


3K 


■EWTJRn 


NIZATION  NAME  and  address 


LuTech,  Inc 
P.0.  Box  1263 
Berkeley,  CA  94701 


".  CONTROLLING  OFFICE  NAME  AND  ADORESS 

Air  Force  Weapons  Laboratory  (NTYE) 
Kirtland  Air  Force  Base,  NM  87117 


10.  PROGRAM  element.  PPCJCCT.  task 
area  A  work  unit  numbers 

64747F/T20&6Q6  ^  r  , 


AGES 


1A.  MONITORING  AGENCY  name  ft  ADDRESS.1'//  different  from  Controlling  Office) 


is.  SECURITY  Class.  (o t  tht$  report) 


Unclassified 


ISa  DECLASSIFICATION.  DOWNGRADING 
SCHEDULE 


<6.  Distribution  statement  (at  i/iit  Rtpcn) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  tty  abetted  entered  In  Block  20.  If  different  from  Report; 


»8.  Supplementary  notes 


9  KEY  WORDS  'Continue  on  revere*  tide  //  neceeeery  end  identify  by  block  number) 

Mode  Suppression 
Parallel -Plate  Simulators 
Modal  Filter 


20  ABSTRACT  fContinuo  on  reveree  tide  It  neceeeery end  identity  by  block  number) 


The  subject  of  this  report  is^he  suppression  or  damping  of  non-TEM  modes  in 
the  transmission-line  type  of  E$P  simulators— +his--suppressfoiv>may  be  accom¬ 
plished  by  introducing  a  spatial  modal  filter  (SMF) ..initially  in  the  output 
cnical  section  of  the  simulator.  v4rf  £his  reporty-wg'  dealj^with  the  electro¬ 
magnetic  onsiderations  of  the  SMF  which  is  essentially  decoupled  from  the 
principal  TEM  mode.  _  / 


00  ,  :«Tn  1473. 


UNCLASSIFIED 


$ECuRitv  Classification  of  this  rage  'Whan  Dtit  Enttrte) 


-  4* 


PREFACE 


We  are  thankful  to  Drs.  K.C.  Chen,  M.G.  Harrison, 
and  Mr.  W.  Kehrer  for  their  encouragement  and  interest 
in  this  problem.  For  many  useful  discussions,  thanks 
are  also  due  to  Dr.  F.M.  Tesche. 


i/ii 


CONTENTS 


Section  Page 

Preface  i 

Contents  1 

List  of  Illustrations  2 

I  Introduction  3 

II  Departure  from  the  Ideal  TEM  Behavior  5 

III  Characteristics  of  non-TEM  Modes  14 

IV  EM  Considerations  of  Spatial  Modal  Filters  18 

V  Estimation  of  Sheet  Impedances  21 

VI  Summary  30 

References  32 


1 


LIST  OF  ILLUSTRATIONS 


Figure 

Page 

2.1 

Geometry  of  a  vertically  polarized,  two 
parallel-plate  transmission-line  type  of 

EMP  simulator 

7 

CM 

• 

CM 

Measured  j Bx/constant i  at  the  center  with 
coordinates  (0,0,0)  in  ALECS  facility 

9 

2.3 

Measured  | Dy/constant I  at  the  center  with 
coordinates* (0 , 0 , 0)  in  ALECS  facility 

10 

2.4 

Measured  |Bx/constant i  at  (0,0,0)  in  the 
model  simulator  at  Harvard 

12 

5.1 

Structure  of  the  volumetric  spatial  modal 
filter 

22 

5.2 

A  typical  TEM  equipotential  and  electric 
field  lines  in  the  conical  transmission  line 
[(b/a)  =  0.857  and  (L/b)  =  4.2;  scale  model 
simulator  with  Z^TEM^  *  82  ft  from  top  plate 
to  ground  plane] 

23 

5.3 

Mode  suppression  sheets  viewed  as  plane 
wave  absorbers 

25 

S.l 

Side  view  of  the  simulator  showing  the 
approximate  locations  of  two  possible  spatial 
modal  filters  (SMF) 

31 

2 


1 


I .  INTRODUCTION 

An  important  class  of  EMP  simulators  is  the  parallel- 
plate  transmission-line  type.  This  type  of  an  electro¬ 
magnetic  structure  is  essentially  an  open  waveguide,  between 
the  conductors  of  which  a  transient  pulse  with  a  planar 
wavefront  must  travel.  In  practice,  however,  the  fields 
obtained  in  the  working  volume  can  depart  from  the  ideal 
TEM  behavior  owing  to  the  excitation  and  propagation  of 
non-TEM  modes.  It  has  been  known  for  some  time  [1  to  5] 
that  a  finitely  wide  parallel-plate  transmission  line  can 
support  and  propagate  TE  and  TM  modes  if  they  are  excited 
for  any  reason.  Since  the  chief  object  of  such  an  EMP 
simulator  is  to  produce  an  EMP  environment  appropriate  to 
a  plane  wave  (outside  the  source  region) ,  it  is  desirable 
to  suppress  the  non-TEM  modes  without  disturbing  the  TEM 
wave.  To  successfully  accomplish  this  suppression,  a 
clear  understanding  of  the  characteristics  of  the  TEM  and 
non-TEM  (i.e.,  TE  and  TM)  modes  is  essential.  Much 
detailed  work  has  been  done  concerning  the  TEM  properties 
of  both  the  parallel-plate  (6  to  9]  and  the  conical 
[10,  11]  transmission  lines.  References  [12,13]  have 
considered  the  conical  transmission  line  as  a  launcher 
and  receptor  of  waves  on  the  cylindrical  transmission  line 
by  introducing  the  concept  of  dispersion  distances,  or 
equivalently,  dispersion  times.  Formuals  have  also  been 
developed  [13]  for  the  TEM  mode  coefficient  in  terms  of 
the  cross-section  fields.  Detailed  calculations  of  TE 
and  TM  modes  of  propagation  are  currently  available  for 
the  two  limiting  cases  of  narrow  (separation  >>  width) 

[1]  and  wide  (separation  <<  width)  [4,5]  plates.  A  para¬ 
metric  study  is  presently  in  progress  [14]  for  the  general 
case  where  the  separation  to  width  ratio  is  not  restricted. 
This  study  will  consider  geometries  of  the  existing  (ALECS 
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and  ARES)  and  future  (ATLAS  I  and  II)  transmission-line 
type  of  EMP  simulator  facilities  at  the  Air  Force  Weapons 
Laboratory,  as  well  as  the  laboratory  model  simulator  at 
Harvard  University.  Comparisons  of  the  electromagnetic 
field  calculations  from  this  study  with  the  available 
experimentally  measured  fields  in  ALECS  and  Harvard's 
model  simulator  is  expected  to  lead  to  an  identification 
of  the  effects  of  the  higher  order  modes.  In  any  case, 
there  exists  a  need  for  suppressing  or  damping  the  non- 
TEM  modes . 

In  Section  II,  the  departure  of  the  measured 
fields  in  the  working  volume  from  the  ideal  TEM  behavior 
will  be  considered,  and  in  Section  III,  available  compu¬ 
tations  of  TE  and  TM  modes  in  the  parallel-plate  region 
will  be  reviewed.  Sections  IV  and  V  deal  with  various 
aspects  of  the  spatial  modal  filter,  and  design  formulas 
are  developed.  The  note  is  concluded  with  a  summary  in 
Section  VI . 
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II.  DEPARTURE  FROM  THE  IDEAL  TEM  BEHAVIOR 

This  class  of  EMP  simulators  operating  in  a  pulsed 
mode  is  a  complex  electromagnetic  structure  to  analyze. 

The  complexity  is  partially  due  to  the  fact  that  the 
input  pulse  contains  a  wide  range  of  frequencies  and, 
consequently,  the  relevant  dimensions  of  the  structure 
ranges  from  a  small  fraction  of  a  wavelength  to  many 
wavelengths.  Because  of  this,  the  simulator  while  operating 
in  a  pulsed  mode  is  a  transmission  line,  a  radiator  and 
an  optical  diffracting  structure,  all  for  the  same  pulse. 

In  reference  [15] ,  the  electromagnetic  characteristics 
of  the  simulator  were  qualitatively  discussed  by  cate¬ 
gorizing  the  frequency  range  of  interest  into  a)  low 
frequencies,  b)  high-frequency  asymptotics,  and  c)  inter¬ 
mediate  frequencies. 

At  low  frequencies,  one  has  near-ideal  conditions 
in  terms  of  simualtion  because  quasi-static  considerations 
apply  and  the  TEM  mode  of  propagation  is  dominant.  The 
main  problem  here  is  to  minimize  the  impedance  dis¬ 
continuity  and  TEM  field  discontinuity  across  the  junction 
between  the  cylindrical  and  the  conical  transmission  lines. 
Such  a  matching  of  the  TEM  modes  at  the  input  and  output 
"bends"  is  achieved  by  reducing  the  dispersion  distance 
[12,13]  . 

At  higher  frequencies,  the  relevant  dimensions  of 
the  simulator,  e.g.,  width  and  spacing  of  the  plate  become 
several  wavelengths  long,  and  ray-optical  considerations 
apply.  Several  canonical  problems  have  been  defined  and 
solved  [16-19]  ,  which  are  useful  in  estimating  the  early- 
time  pulsed  fields  associated  with  conducting  wedges 
and  thus  reduce  waveform  distortion.  Some  experimentally 
measured  data  [20]  is  available  concerning  the  early-time 
fields  indicating  the  (1/r)  variation  of  the  spherical 
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wave  launched  by  the  conical  line  followed  by  edge  diffrac¬ 
tion  and  specular  reflections.  Some  of  the  problems  here 
are  large  dispersion  times,  edge  diffractions  and  ex¬ 
cessive  ripple  in  the  measured  data. 

However,  it  is  the  intermediate  frequencies  that 
present  the  most  serious  problems.  Relevant  distances 
now  become  comparable  to  wavelengths  rendering  both  the 
quasi-static  and  the  ray-optic  considerations  inapplicable. 
In  this  note,  we  focus  our  attention  on  the  behavior  of 
the  simulated  electromagnetic  fields  in  this  frequency 
regime  by  taking  a  superposition  of  the  TEM  and  non-TEM 
modes.  The  non-TEM  modes  are  those  supported  in  a 
finitely  wide  two-parallel-plate  open  waveguiding 
structure,  which  are  to  be  contrasted  with  the  familiar 
propagating  modes  of  a  closed  rectangular  waveguide  or 
a  closed  circular  coaxial  transmission  line.  By  taking 
this  view,  one  can  define  the  problem  to  be  that  of 
damping  or  suppressing  the  non-TEM  modes  without  sig¬ 
nificantly  disturbing  the  desired  TEM  mode  of  propagation. 

A  typical  geometry  of  this  class  of  EMP  simulator 
is  shown  in  Figure  2.1.  This  figure  shows  a  vertically 
polarized  parallel-plate  transmission-line  type  of 
simulator.  By  virtue  of  symmetry  and  practicality ,  it 
suffices  to  construct  the  symmetric  half  of  the  structure 
above  the  horizontal  symmetry  plane  where  a  wide  conducting 
plane  is  placed.  Accordingly,  Figure  2.1  shows  the  side 
view  of  the  simulator  comprising  the  top  plate  of  width  2a 
at  a  height  b  above  the  ground  plane.  The  plate-to-plate 
separation  (for  a  symmetric  situation  with  two  identical 
parallel  plates)  is  then  2b  and  the  lengths  of  the  conical 
and  the  parallel-plate  transmission  lines  are  designated  L 
and  w  respectively.  With  reference  to  this  figure,  the 
nonzero  components  of  electric  and  magnetic  fields  for  the 
various  modes  of  propagation  in  the  parallel-plate  trans¬ 
mission-line  region  are: 
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plate  transmission-line  type  of  EMP  simulator 


a)  TEM  mode 


£  ,  E  ,  H 

x  y  x 

,  H 
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with 

Ez  =  Hz  =  ° 
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V  *x  '  *y  ' 
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,  H 
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with 
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TE 
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with 

e2  =  0 

All 

of 

the  above  field  quantities  are 

functions  of  position 

in 

the 

transverse  plane 

and 

frequency . 

We  shall  now 

summarize  certain  measurements  made  in  the  ALECS  facility 
that  are  relevant  in  terms  of  identifying  the  departure 
from  the  ideal  TEM  behavior.  CW  measurements  [21]  made 
in  the  working  volume  of  the  ALECS  facility  have  detected 
what  has  been  referred  to  as  the  "notch  problem". 

Specifically,  the  transfer  function  from  the  input  voltage 
to  the  measured  fields,  when  appropriately  normalized 
and  plotted  as  a  function  of  frequency,  displays  signifi¬ 
cant  notches  at  certain  frequencies.  For  example,  Figure  2.2 
shows  the  magnitude  of  the  normalized  transverse  magnetic 
field  H  measured  at  the  center  (0,0,0)  and  the  notch  is 

X 

seen  to  appear  at  ~25  MHz,  with  roughly  a  ±30%  ripple  at 
higher  frequencies.  Figure  2.3  is  the  normalized  principal 
electric  field  E  measured  at  the  same  location.  It  is 
noted  that  the  principal  electric  field  does  not  display  a 
notch  behavior  at  the  same  location,  but  perhaps  a  small 
enhancement.  Correspondingly,  the  experimental  measurements 
in  the  scale  model  simulator  at  Harvard  University  have  repro¬ 
duced  the  notch  behavior  in  the  field  quantities.  The  descrip¬ 
tion  of  the  experimental  setup  along  with  the  measured  fields 
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Figure  2.3  Measured  |Dy/constant|  at  the  center  with 
coordinates  (0,0,0)  in  ALECS  facility 

(This  figure  is  reproduced  here  from  reference  M> 


in  the  model  simulator  are  well  documented  in  references 
[22]  and  [23]  .  A  preliminary  measurement  of  the  trans¬ 
verse  magnetic  field  as  a  function  of  frequency  at  a 
fixed  location  of  (0,0,0)  was  provided  by  Blejer  [24] 
and  is  shown  in  Figure  2.4.  This  measurement  also 
displays  a  sharp  null  in  the  magnetic  field  at  -264  MHz 
at  the  measurement  location.  The  null  is  caused  by  a 
cancellation  of  the  transverse  magnetic  field  of  the 
principal  TEM  mode  by  that  of  the  first  higher  TM  mode, 
and  the  frequency  of  264  MHz  on  the  model  simulator 
corresponds  very  closely  to  the  -25  MHz  notch  observed 
in  the  ALECS  facility.  Furthermore,  in  both  configurations 
(ALECS  and  the  model  simulator) ,  the  frequencies  where 
one  would  expect  a  higher  order  mode  to  be  excited 
corresponds  approximately  to  the  relevant  dispersion 
distance  "d"  becoming  equal  to  a  half  wavelength. 
Considering  a  direct  path  from  source  apex  to  load  apex 
and  another  signal  path  along  the  edge  of  the  top  plate, 
d  is  given  by 


d 


2  2 
a  +  L 


(2.1) 


where  X  is  the  wavelength  corresponding  to  a  frequency 
where  one  may  expect  a  higher  order  TM  mode  excitation, 
and  the  other  variables  in  Eq.  (2.1)  are  illustrated  in 
Figure  2.1.  Substituting  the  values  for  the  various 
dimensions,  we  find  that  the  frequency  of  expected  exci¬ 
tation  of  higher  TM  mode  is  ~25  MHz  for  ALECS  and 
-256  MHz  for  the  model  simulator  at  Harvard.  These 
compare  well  with  the  experimentally  observed  values  of 
-24  MHz  for  ALECS  and  ~264  MHz  for  the  model  simulator. 

It  can  be  argued  that  the  notch  is  due  to  the  superposition  of 
the  TEM  mode  and  higher  order  TM  and  TE  modes.  The  non-TEM 
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modes  being  evanescent  at  either  end  of  the  simulator,  as 
expected,  all  of  the  input  impedance  measurements  are  domi¬ 
nated  by  the  principal  TEM  mode.  Based  on  these  observa¬ 
tions,  one  may  conclude  that  the  notch  occurs  at  the  center 
(0,0,0)  due  to  the  symmetry  in  the  simulator  as  a  result 
of  near  cancellation  of  the  magnetic  field  of  the  TEM  mode 
with  its  counterpart  in  the  first  TM  mode.  We  do  not 
expect  a  near  cancellation  in  the  electric  field  at  the  same 
location  for  the  following  reasons.  The  dominant  TEM  mode 
is  a  traveling  base  because  of  the  notched  load  at  the  end 
whereas  the  TE  and  TM  modes  are  standing  waves  owing  to 
reflections.  The  magnetic  field  (H^  component)  in  the  TM 
modes  propagating  in  the  positive  and  negative  2  directions 
add  and  their  sum  cancels  the  Hx  of  the  principal  TEM  mode 
resulting  in  the  notch  behavior.  However  the  electric  field 
in  the  TM  modes  propagating  in  the  positive  and  negative  2 
directions  cancel  from  each  other  and  thus  one  does  not 
expect  a  notch  at  (0,0,0)  in  the  electric  field.  This  is 
consistent  with  the  experimentally  measured  electric  field 
(Ey)  data  shown  in  Figure  2.3.  Furthermore,  near  the  top 
plate  both  and  E^  in  the  TM  mode  will  have  the 

reversed  sign  from  their  values  at  the  ground  plane,  leading 
to  no  cancellation  in  either  of  these  field  quantities  with 
their  counterparts  in  the  TEM  mode.  Vie  also  note  that  the 
frequency  where  the  notch  occurs  is  also  predictable  from 
the  dispersion  times/distances  calculations. 


13 


III.  CHARACTERISTICS  OF  NON-TEM  MODES 

An  open  waveguide  formed  by  two  perfectly  con¬ 
ducting  parallel  plates  of  finite  width  can  suDoort,  in 
addition  to  the  dominant  TEM  mode,  higher-order  TE  and 
TM  modes.  These  higher  modes  have  suitably  complex 
propagation  constants  which  account  for  the  radiation  loss 
of  the  modes.  Because  of  radiation,  power  flow  and 
stored  energy  are  not  confined  to  the  inner  region  of 
the  structure.  As  the  wave  (a  superposition  of  modal 
distributions)  propagates,  there  is  a  continuous  "leaking" 
of  energy  from  the  open  waveguiding  structure. 

In  closed  waveguides,  the  higher  modes  are  the 

solutions  of  source-free  time  harmonic  Maxwell's  equations 

characterized  by  axial  field  variation  of  the  form 

exp  [  j  (oit-z-k  )  ]  .  Each  one  of  these  modes  with  real 

propagation  constants  (k  )  below  cutoff  satisfies  all 

z 

of  the  boundary  consitions.  Above  cutoff  k2  is  pure 
imaginary  and  the  modes  are  evanescent.  Also,  these 
modes  possess  finite  energies  in  a  cross  section  and 
display  orthogonality  and  completeness  properties. 

In  contrast,  the  modes  on  an  open  waveguide  form 
a  discrete  spectrum  of  leaky  modes  satisfying  Maxwell's 
equations  and  boundary  conditions.  However,  any  cross 
section  extends  to  infinity  and  the  fields  can  grow 
without  bounds  at  large  distances  from  the  simulator  cross 
section.  These  leaky  modes  on  open  waveguides  do  not 
form  a  complete  orthogonal  set  and  have  to  be  supplemented 
by  a  continuous  spectrum.  The  variation  of  an  arbitrary 
mode  can  be  separated  into  transverse  and  longitudinal 
parts  as  follows: 

**  ~  —  £  Z 

F(x,y,z,s)  =  f(x,y,s)  e  (3.1) 
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where  F  denotes  an  arbitrary  field  component,  and  f 
denotes  variation  in  the  transverse  (x-y)  plane.  The 
exponential  factor  displays  axial  (z  coordinate)  variation 
with  5  (Laplace  transformation  variable  corresponding 
to  the  z  coordinate)  being  the  longitudinal  complex 
wavenumber.  The  complex  frequency  is  denoted  by  s,  and 
it  is  the  Laplace  transformation  variable  corresponding 
to  the  time  t  .  Setting, 

2  ,2,2.  2  ,2  2,  ,,  ,, 

p  =  (s  /c  )  -  ?  =  (y  -  c  )  (3.2) 


we  have. 


(3.3) 


To  exhibit  the  cutoff  behavior,  consider  s  =  jw  , 
y  =  jk  and  5  -  jkz 

S  =  jkz  =  t>/( jk)2  +  (±jp)2  =  ±j  >/k2  +  p2  (3.4) 

or,  kz  =  ±  Vk2  +  p2  -  (3  +  ja)  (3.5) 

In  the  integral  equation  formulation  [1] ,  complex  p  , 
multiplied  by  a  factor  with  linear  dimension  appears  in 
the  argument  of  special  functions.  In  order  to  compute 
the  various  field  quantities,  one  looks  for  the  roots  of 
transcendental  complex  functions  of  p  in  the 
p-plane.  These  functions  are  obtainable  in  closed  form 
for  the  limiting  cases  of  narrow  and  wide  plates.  Denoting 
a  typical  root  by  pn  m  ,  the  corresponding  longitudinal 
complex  wavenumber  for  this  arbitrary  mode  becomes 


+ 


~~1 

^n,m 


8  +  j 

n,m  J 


ot 


n  ,m 


(3.6) 
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Pn/In  values  of  interest  are  approximately  pure  imaginary 
with  relatively  small  negative  real  parts.  Let 


:>  „  =  u  +  j  v  ;  (  !u  |  <<  'v  |) 

n,m  n,m  J  n,m  1  n,!!!1  1  n,m'  ' 


so  that 


k  =  ±^l(k2  +  u2  -v2)  +  2j  u  v  I  (3.7) 
z  _  *  n,m  n,m  J  n,m  n,m 

n  f  in 


If  p  was  pure  imaginary  (i.e.,  u  =0)  ,  the 

n,m  n,m 

cutoff  behavior  depends  upon  the  value  of  the  factor 

2  2  2 
±  V k  -  v  .  Non-TEM  modes  do  not  propagate  when  v 


n,m 
2 

exceeds  k  which  makes  k  pure  negative  imaginary 


n,m 


v2  >  k2 
n,m 


“j  I* 


and  the  resulting  argument  of  exp(-j  k  z)  negative 

real.  The  cutoff  behavior  for  the  special  case 

when  [s=jaj»  y=jk  and  p  =j  v  ]  can  be  summarized  as 
J  *n,m  J  nm 

follows.  For  these  special  conditions,  the  longitudinal 

wavenumber  k z  is  given  by 

[with  S„  _  >  0  and  a  m  =  0]  for  propagation 
n,m  n,m  n,m  r  c  r 

'n,m  I  8  +ja  [with  „>  0  and  a  m  <  0]  for  evanescence 

v  n,m  n,m  n,m  n,m 

(3.8) 

Correspondingly,  the  longitudinal  variation  of  the  fields 
exp(~jkZn  ^z)  exhibits  the  transition  from  propagation 
into  evanescence.  However,  the  complex  singularities 
(pn  m)  do  have  small  negative  real  parts,  and  the  transition 
near  cutoff  is  less  abrupt.  In  the  regime  of  evanescence,  the 
longitudinal  wavenumber  kZn  m  becomes  largely  pure  negative 
imaginary  with  a  small  real  part.  Consequently,  the 


exponential  longitudinal  variation  exp(-jk. 


z)  exhibits 


'n  ,m 


evanescence . 

Furthermore,  a  waveguiding  structure  formed  by  two 
finitely  wide  plates  has  two  planes  (x=0  and  y=0)  about 
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which  the  fields  can  be  symmetric  or  antisymmetric  [25] . 
Such  symmetry  decomposition  can  be  represented  symbolically 
by  indexing  the  p-plane  singularities  as 

(±,±,E  or  H) 
pn,m 


where  n,m  correspond  to  the  field  variations  in  x 
and  y  directions.  The  two  ±  signs  in  the  superscript 
indicate  symmetry  or  antisymmetry  along  the  transverse 
x  and  y  directions,  and  E  or  H  denote  TM  or  TE 
mode.  Thus,  such  a  representation  of  the  p-plane  singu¬ 
larities  uniquely  specifies  a  mode  and  its  symmetry 
properties.  Consistent  with  the  above  notation,  the 


modes  themselves  can  be  denoted  by  E 


{±rt)  (or  TM(±'±)) 
n,m  n,m 


/  \  /X  ,  XXX  11  ,  XXX 

and  HU'i)  (or  TE„'±)  ). 
n  ,m  n,m 

These  modes  in  open  waveguides  formed  by  a  pair  of 
finitely  wide  parallel  plates  are  found  by  formulating 
two  different  scalar  integral  equations  (of  the  first 
kind)  for  the  current  and  charge  densities  on  the 
plate.  Under  certain  approximations  (separation  >>  width 
or  vice  versa) ,  the  integral  equations  can  be  solved 
analytically  by  first  tranforming  them  into  a  Fredholm 
integral  equation  of  the  second  kind  [1,4,5]  and  using 
perturbation  techniques. 

It  is  not  our  intention  to  show  detailed  modal 


distributions  of  these  higher  modes  in  this  section  but 
to  point  out  the  methods  employed  and  availability  of 
field  plots.  Specific  calculations  of  the  TM  modes  in  a 
geometry  corresponding  to  the  model  EMP  simulator  at 
Harvard  are  currently  in  progress  for  later  comparisons 
with  the  measured  data.  Such  comparisons  are  expected  to 
lead  to  unambiguous  identification  of  modes. 
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IV.  ELECTROMAGNETIC  CONSIDERATIONS  OF  SPATIAL  MODAL  FILTERS 

The  object  of  introducing  a  spatial  modal  filter 
is  to  load  or  damp  the  non-TEM  modes  (i.e.,  TE,  TM) 
without  significantly  disturbing  the  TEM  modes.  One 
method  that  has  been  experimentally  implemented 
in  the  context  of  TEM  cells  [26]  with  some  degree  of 
success  is  by  inserting  RF-absorbing  material  [27,28] 
along  the  conducting  walls  of  the  waveguiding  structure. 

This  technique  lacks  an  analytical  basis  and  also  does 
not  fully  exploit  the  uniform  characteristics  of  the 
dominant  TEM  mode  for  damping  the  non-TEM  modes.  If  we 
recognize  the  fact  that  all  propagating  modes  have  spatial 
properties,  one  can,  hence,  take  advantage  of  certain 
spatial  properties  to  selectively  load  the  modes  [29] . 

In  a  two-parallel-plate  transmission-line  type  of  EMP 
simulator,  all  of  the  important  characteristics  of  the 
principal  TEM  mode,  both  in  the  parallel-plate  region  [9] 
and  the  conical  plate  region  [10,11],  are  well  known  and 
documented.  The  quantities  of  interest  are  the  charac¬ 
teristic  impedance,  equipotential  and  field  distributions, 
and  field  uniformity.  Later  in  this  section  we  shall  use 
a  knowledge  of  these  quantities  in  developing  a  spatial 
modal  filter.  Such  filters  can  take  various  forms  and  we 
dimensionally  categorize  them  below: 

A.  Zero  Dimensional 

A  zero  dimensional  spatial  modal  filter  will 
essentially  act  like  a  directional  coupler  in  a  transmission 
line  that  couples  to  backward  propagating  modes.  In  the 
case  of  a  simulator  under  consideration,  this  type  of  SMF 
can  take  the  form  of  two  receiving  or  parasitic  antennas 
(equivalent  dipoles,  loops)  suitably  oriented  on  the  top 
plate  or  the  ground  plane  so  that  the  TEM  mode  is  uncoupled. 
The  sensitivity  of  damping  of  the  non-TEM  modes  is  governed 
by  local  field  ratios. 
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B .  One  Dimensional 

One-dimensional  spatial  modal  filters  can  be  imple¬ 
mented  by  having  longitudinal  slots  both  in  the  top  plate 
and  the  ground  plane.  The  longitudinal  slots  (directed  along 
the  2-axis  of  Figure  2.1  in  the  parallel-plate  region  and 
directed  along  the  spherical  radial  coordinate  in  the  conical 
region)  will  be  resistively  loaded  by  transverse  resistors 
(x-direction  in  Figure  2.1).  By  definition,  the  only  mode 
with  longitudinal  component  of  magnetic  field  (Hz)  is  the 
TE  mode  and,  hence,  this  mode  has  principally  transverse 
currents  (J^)  which  couple  to  the  transverse  resistors. 
Consequently,  this  form  of  SMF  is  essentially  uncoupled 
from  the  principal  TEM  or  higher  TM  modes  while  loading 
the  higher  TE  modes. 

C.  Two  Dimensional 

An  extension  of  the  one-dimensional  SMF  is  to 
introduce  "loading  sheets"  that  are  comprised  of  a  two- 
dimensional  (transverse  and  longitudinal)  array  of 
resistors.  In  the  top  plate  and  the  ground  plane,  one 
would  have  only  transverse  resistors,  but  in  the  space 
away  from  the  simulator  plates,  the  loading  sheet  will 
be  two  dimensional.  The  sheet  will  be  located  along  an 
equipotential  surface  of  the  principal  TEM  mode.  This 
implies  that  the  sheet  is  electromagnetically  invisible 
to  the  principal  TEM  mode  while  coupling  and,  hence, 
damping  the  non-TEM  modes. 

D.  Three  Dimensional 

The  two-dimensional  loading  sheet  described  above 
can  be  repeated  to  fill  some  portion  of  the  volume  between 
the  top  plate  and  the  ground  plane  resulting  in  a  volumetric 
suppressor  of  non-TEM  modes. 

There  are  certain  constraints  to  be  placed  on  the 
choice  of  the  volume  of  space  wherein  such  a  filter  can  be 
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placed.  For  example,  the  filter  should  not  be  close  to  the 
working  volume  of  the  simulator  to  avoid  any  coupling  to 
the  test  object.  Furthermore,  the  higher-order  TE  and  TM 
modes  become  evanescent  at  a  certain  distance  away  from 
the  terminator  in  the  output  section.  For  this  reason, 
little  is  gained  by  placing  the  volumetric  SMF  near  the 
termination.  Such  considerations  indicate  that  the  SMF 
should  be  placed  under  the  output  bend  extending  toward, 
but  not  close  to,  the  terminator. 

In  Section  V,  we  develop  an  approximate  analytical 
basis  for  estimating  the  sheet  resistors. 
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V.  ESTIMATION  OF  SHEET  IMPEDANCES 

In  this  section,  we  derive  relations  useful  in 
estimating  the  longitudinal  and  tranverse  resistor  values 
in  the  individual  sheets  of  a  volumetric  suppressor.  A 
cross-sectional  view  of  the  SMF  is  sketched  in  Figure  5.1. 
Essentially  the  mode  filter  would  consist  of  several 
layers  of  resistor  arrangement,  wherein  each  layer 
consists  of  a  two-dimensional  array  of  resistors  along 
the  axial  and  tranverse  directions.  Each  layer  would 
coincide  with  a  TEM  equipotential  surface  in  a  conical 
transmission  line.  The  equipotential  and  field  calcula¬ 
tions  in  a  conical  transmission  line  have  been  reported 
[10,11],  and  they  are  based  on  stereographically  pro¬ 
jecting  the  two-conical-plate  line  into  a  cylindrical 
line  of  two  circular  arcs  on  two  different  circles.  The 
curved-cylindrical-plate  problem  is  later  solved  by  the 
method  of  conformal  mapping  for  the  TEM  quantities. 

The  physical  quantity  that  is  relevant  for  the  SMF 
design  is  the  TEM  equipotential  surface,  an  example  of 
which  is  shown  plotted  in  Figure  5.2.  In  this  figure, 
the  top  plate  is  held  at  a  potential  of  VQ  with  the 
ground  plane  as  the  reference.  The  equipotentials  (magnetic 

field  lines)  are  shown  in  steps  of  0.1  V  and  the  stream 

o 

lines  (electric  field  lines)  in  steos  of  0.1  U  .  U  is  the 

o  o 

total  current  flowing  on  half  of  the  top  plate.  The 

individual  resistive  sheets  of  the  volumetric  suppressor 

will  be  made  to  coincide  with  the  calculated  surfaces. 

In  estimating  the  values  of  the  transverse  (Rfc) 

and  the  longitudinal  ( )  resistors,  it  is  useful  to 

view  the  two-dimensional  array  of  resistors  on  any  given 

equipotential  surface  as  a  resistive  sheet  with  an  impedance 

of  Z  .  We  now  consider  the  following  two  canonical 
s 

problems  useful  in  estimating  R^  and  Rfc  .  The  two 
problems  are  posed  and  analyzed  by  viewing  the  mode 
suppressor  sheets  as  plane  wave  absorbers. 
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Cross  Sectional  View 


volumetric  spatial  modal  filter 
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Figure  5.2  A  typical  TEM  equipotential  and  electric 

field  lines  in  the  conical  transmission  line 
[ (b/a)  =  0.357  and  (L/b)  =  4.2;  scale  model 
simulator  with  zjTEM)=  8 2f2  from  top  plate 
to  ground  plane] 
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A.  TM  Waves  Incident  on  the  Resistive  Sheet 

Consider  a  TM  wave  incident  on  the  resistive  sheet 
( Z  in  f2)  ,  as  shown  in  Figure  5.3a.  Denoting  the  induced 
surface  current  density  on  the  sheet  by  Jg(A/m),  the 
boundary  conditions  on  the  tangential  components  of  the 
electric  and  magnetic  field  lead  to 

E.__  =  E;  sin(a)  -  E  sin(a)  =  E.  sin  (a)  (5.1) 

can  1  r  c 

The  difference  in  the  tangential  magnetic  field  is 

tan  top  bottom  l  r  t  s  tan  s 
J  is  positive  and  is  directed  from  left  to  right  in 
Figure  5.3a.  Note  that  the  free  space  characteristic 
impedance  Zq  is  given  approximately  by 


Using  Eq.  (5.3)  in  Eq.  (5.2),  and  then  making  use  of 
Er  =  E^  -  Efc  from  Eq.  (5.1) ,  we  can  get 

Ei-  2 

r~  ...  _ s _  ,c  4  \ 

E”  "  2  7,  +  sin  (a) 

1  so 

For  fixed  values  of  E^  and  a  ,  the  power  dissipated  in 
the  sheet  per  unit  area. 


P  = 


E.  J 
tan  s 


-  <ELn/ZS>  *  <E? 

E2  sin2 (a)  (  2  Zg  ^ 

=  z;  [. 2  Zs  +  ZQ  sin  (a)  y 


(5.5) 


Figure  5.3  Mode  suppressor  sheets  viewed  as  plane  wave  absorbers 
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Now  the  power  dissipated  as  heat  from  the  incident  wave 
per  unit  area  normal  to  the  incident  wave  is 

■2  i 

(5.6) 


Pw  =  P/sin 


in(«,  -  f) 


where  f  is  the  fractional  power  lost,  and  is  given  by 

Z  Z  /  2  Z  \2 

f  _2  sin  (a)  Ze  Sln(at)  l  T~Z  +~1  sin '(a)  ) 

E.  s  \  s  o  / 


(5.7) 


Setting 


n  = 


Z  sin  (a) 
o 


(5.8) 


we  have 


2n 


1  +  2n 


(5.9) 


We  may  now  maximize  the  fractional  power  lost  by  requiring 

af  "  0  ?  leading  to  nopt  =  1/2 

This  gives  an  optimum  value  for  the  sheet  impedance  of 

(5.10) 


Z  sin  (a) 
2(opt)  __  o 


s 


Equation  (5.10)  is  now  useful  in  estimating  the  longitudinal 
resistors  which  carry  the  current  for  the  TM  mode  of 

propagation. 
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B.  TE  Waves  Incident  on  the  Resistive  Sheet 


We  can  use  an  entirely  similar  analysis  as  in  the 
TM  case  for  the  TE  case.  Consider  a  TE  wave  incident  on 
a  resistive  sheet  of  impedance  Zg  in  Q  at  an  angle  a  as 
shown  in  Figure  5.3b.  Let  the  induced  surface  current  on 
the  resistive  sheet  be  denoted  by  Jg(A/m).  The  boundary 
conditions  on  the  tangential  fields  are 


tan 


=  -E.  -  E  =  -E. 


(5.11) 


The  difference  in  the  tangential  magnetic  field  is 


AHtan  =  Htop  “  Hbottom  =  ~Hi  sin(a)  +  Hr  sin(a)  +  Ht  sin  (a) 


top 
=  J  =  (E 


tan/V. 


(5.12) 


Jg  is  positive  and  directed  into  the  plane  of  the  paper  in 
Figure  5.3b.  As  before,  using  Eq.  (5.3)  in  (5.12)  and  later 
using  Eq.  (5.11),  we  can  get 


E.  2  Z  sin(a) 

_t  -  s _  (5 

E.  2  Z  sin (a)  +  Zo 
is  o 

Once  again,  the  power  dissipated  in  the  sheet  per  unit 
area  is 


P  =  E ,  J 

tan  s 


(EL/Zs>  '  (El/Zs> 


■S( 


Z  sin(a) 


2  Z„  sin (a) 
s 


+  z. 


(5.14) 


Using  Eq.  (5.6),  the  fractional  power  lost  is  now  given 
by 


Zo  P  Zo  /  2  Zs  sin(a)  \2 

f  =  ^2  sin (a)  =  Zg  sin(a)  ^2  Zg  sin(a)  +  Z ~J  (5,15) 

-h  [rr^]  <5-16’ 

where  v  =  (Zg  sin(a)/ZQ).  We  may  now  maximize  the 
dissipated  power  by  requiring 


df 

dv 


0 


leading  to  vQpt  =  1/2 


(5.17) 


This  gives  an  optimum  value  for  the  sheet  impedance  of 


2 (opt)  ..  2o 
s  "2  sin(ct) 


(5.18) 


Once  again,  we  note  that  for  the  TE  case,  only  the  trans¬ 
verse  resistors  carry  current  and,  hence,  Eq.  (5.18) 
applies  to  the  transverse  resistors  Rfc  . 

In  either  case  of  TM  or  TE  waves  incident  on  the 
resistive  sheets,  the  values  of  the  resistors  R^  and 
Rfc  are  dependent  on  the  angle  of  incidence  of  these  waves 
onto  the  sheets.  In  a  simulator  configuration,  this  angle 
(a)  is,  of  course,  a  variable  quantity  and,  consequently, 
an  experimental  optimization  of  the  particular  values  of 
the  resistors  is  inevitable. 

However,  at  high  frequencies,  where  ray-optic 
considerations  apply,  one  can  estimate  the  angle  a  by 
considering  a  typical  ray  path.  It  is  estimated  that  the 
present  high  frequency  ripple  of  +30%  can  be  reduced  to 
within  +10%.  At  intermediate  frequencies,  the  angles  are 
harder  to  estimate,  but  effective  removal  of  energy  from 
the  non-TEM  modes  results  from  multiple  passes.  It  is 
likely  that  for  E  modes  (TM  case) ,  the  ratio  of  Ez  to 
E  or  e/  may  define  an  effective  angle  a  .  A  future 

x  /  y 

memo  will  address  these  issues  specifically. 
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For  the  experimental  evaluation,  initially,  a 
typical  angle  may  be  chosen  in  designing  the  filter  and, 
later,  experimentation  around  these  values  will  determine 
the  final  values. 

In  successfully  implementing  these  concepts  in  a 
given  experimental  situation,  there  will  be  several 
associated  problems  and  considerations  unique  to  the 
experiment.  For  instance,  in  a  model  simulator  experiment, 
the  power  levels  are  low  and,  hence,  energy  dissipation 
in  the  resistors  is  not  a  serious  consideration.  However, 
on  a  full-scale  facilitiy  (e.g.,  ALECS),  it  may  become 
necessary  to  place  the  resistors  in  plastic  tubes  filled 
with  oil  (if  required) .  Also,  in  an  on-site  configuration, 
the  mechanical  problems  of  supporting  and  holding  the 
sheet  along  equipotential  contours  are  much  more  severe. 
There  are  questions  like  how  far  outward  in  the  cross 
section  should  the  sheets  extend.  Such  considerations 
and  a  detailed  design,  fabrication  and  evaluation  of  the 
proposed  volumetric  suppressor  will  form  the  subject  of 
a  separate  report. 
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VII.  SUMMARY 

A  two-parallel-plate  transmission-line  type  of 
EMP  simulator  has  been  historically  employed  to  propagate 
an  electromagnetic  pulse  with  a  plane  wavefront.  Since 
the  pulse  consists  of  a  wide  frequency  spectrum,  the 
structures  that  have  been  built  propagate  modes  other  than 
the  principal  TEM  mode.  Consequently,  techniques  of 
non-TEM  mode  suppression  is  an  important  element  in  the 
advancement  of  simulator  technology.  In  this  note,  we 
address  the  problem  of  mode  suppression  and,  in  particular, 
the  electromagnetic  considerations  of  a  volumetric 
suppressor..  This  suppressor  is  the  result  of  one  of 
several  techniques  of  mode  suppression  and  believed  to 
be  the  most  efficient.  Empirical  relations  are  developed 
for  computing  the  resistor  values  required  in  the  fabri¬ 
cation  by  viewing  the  suppressor  as  a  plane  wave  absorber. 

An  important  feature  of  the  SMF  is  its  location. 
Figure  6.1  shows  the  side  view  of  a  parallel-plate  EMP 
simulator,  indicating  the  approximate  locations  of  two 
possible  SMF.  The  SMFs  have  to  be  located  away  from  the 
two  apexes  where  the  non-TEM  modes  are  evanescent.  It  is 
also  essential  to  place  the  SMF  so  as  to  minimize  any 
coupling  effects  with  the  test  object. 

Finally,  although  our  present  interest  is  in  EMP 
simulators,  it  is  recognized  that  the  mode  suppression 
considerations  presented  here  are  fairly  ger.  .al  and  can 
apply  to  other  types  of  multimoded  structures  (e.g.,  TEM 
cell)  where  the  non-TEM  modes  are  also  undesirable,  and 
due  to  the  lack  of  radiation  losses  can  cause  even  more 
severe  problems. 
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5 


initial  SMF  in  the  output 
section 


possible  second  SMF  in 


Figure  6.1  Side  view  of  the  simulotor  showing  the  approximate 
locations  of  two  possible  spatial  modal  filters  ( SMF) 
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